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Molecular electronic devices currently serve as a platform for studying a variety of physical 
phenomena only accessible at the nanometer scale. One such phenomenon is the highly correlated 
electronic state responsible for the Kondo effect, manifested here as a "Kondo resonance" in the 
conductance. Because the Kondo effect results from strong electron-electron interactions, it is not 
captured by the usual quantum chemistry approaches traditionally applied to understand chemical 
electron transfer. In this review we will discuss the origins and phenomenology of Kondo resonances 
observed in single molecule devices, focusing primarily on the spin-i Kondo state arising from a 
single unpaired electron. We explore the rich physical system of a single-molecule device, which 
offers a unique spectroscopic tool for investigating the interplay of emergent Kondo behavior and 
such properties as molecular orbital transitions and vibrational modes. We will additionally address 
more exotic systems, such as higher spin states in the Kondo regime, and we will review recent 
experimental advances in the ability to manipulate and exert control over these nanoscale devices. 
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Research in the field of molecular electronics has been 
driven by an assortment of promising technological ad- 
vancements, both fundamental and applied in nature. 
The discrete building blocks used to construct modern 
electronic components will ultimately limit the extent 
to which these devices can be scaled down in size. De- 
vice miniaturization may be taken to new extremes by 
fabricating components that utilize a single molecule or 
molecular complex as the active element. The realization 
of such devices has presented a challenge to our ability to 
manipulate and exert control over these nanoscale struc- 
tures, and has provided new insights and an expanded 
understanding of the way in which a molecule may inter- 
act with solid state components. 

Early theories of electrical conduction through 
molecules were presented in the 1940's by Robert Mul- 
likcn and Albert Szcnt-Gyorgi. However, a practical 
foundation for the premise of a device operating with 
a single, custom designed molecule is widely credited as 
having begun with the theories of Aviram and Ratner, 
put forth in 1974pi They proposed using a donor-acceptor 
molecule as a molecular diode. Shortly thereafter came 
an upsurge of papers presenting more detailed analyses 
of molecular conduction and designs for molecular wires, 
switches, and logic circuits. The elementary structure of 
all of these devices, regardless of the application, con- 
sists of a molecule or molecular complex attached to two 
metallic electrodes functioning as source and drain con- 
tacts. The incorporation of a third, capacitively cou- 
pled electrode, typically referred to as a gate, allows for 
the formation of a 3-terminal transistor-like device. This 
structure is then analogous to the field effect transistor 
(FET), the basis for nearly all modern electronic devices. 



The experimental realization of even the most basic of 
such designs first required the development of new fabri- 
cation techniques enabling means of contacting individ- 
ual molecules via metallic electrodes. These techniques 
led to early transport studies, primarily directed at inves- 
tigating conduction through molecules via two terminal 
I—V measurements. The field has taken sizable steps for- 
ward recently thanks to collaborations between chemists 
synthesizing tailor-made complexes, and physicists incor- 
porating these complexes into nanoscale devices. 

This review focuses primarily on Kondo physics in 
single-molecule electronic devices. As we explain in 
more detail below, the Kondo effect, first noticed in the 
1930s in bulk metals containing magnetic impurities, is 
a consequence of strong electron-electron interactions 
and correlations. It therefore lies outside the scope of 
single-particle quantum chemistry descriptions, while 
profoundly affecting electron transport. The Kondo 
effect leads to a pronounced zero bias Kondo reso- 
nance in the conduction that would be absent in a 
single-particle picture. We first introduce the fabrication 
techniques required to produce single-molecule electronic 
devices. After an overview of the physics of electron 
transport in such nanodevices, we discuss the origin 
of the Kondo effect. Implementations of the Kondo 
effect in molecules arc presented together with the 
means by which these systems can be probed in order 
to better understand the correlation effects. We also 
discuss the way in which experimental observations 
deviate from commonly used models of the Kondo state. 
Finally, we address scenarios beyond the single channel 
spin-1/2 system and possible future lines of investigation. 
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FABRICATION 

Single molecule devices, serving here as a platform for 
the study of the Kondo effect, are often described as a 
double barrier tunneling structure in which a molecule is 
connected to metallic source and drain electrodes by con- 
tacts that each act as tunneling barriers. This arrange- 
ment can be accomplished by more than one method. 

In a scanning tunneling microscope (STM), the tun- 
neling current is measured between a sharp metallic tip 
(capable of controlled movement in the x, y, and z direc- 
tions) and a conducting surface. This tool enables one 
to map the density of states (DOS) of a surface covered 
with dispersed molecules. When a molecule of interest 
is found, situating the molecule between the sharp tip 
and the conducting substrate establishes a double bar- 
rier junction (Fig. QJi). Electrons can tunnel between 
the substrate (including the molecule) and the atom on 
the STM tip that is situated closest to the substrate sur- 
face. The absence of a gate electrode limits the STM to 
two terminal measurements. Kondo physics in STM ex- 
periments has been observed for more than a decade*^ 
with much recent progress being made in molecular sys- 
tems, as described below. 

Alternately, planar single-molecule devices may be fab- 
ricated using "breakjunction" techniques. Arrays of 
nanowires are lithographically defined on an insulated 
substrate. A dilute solution of the molecules of interest is 
spread on the substrate surface, with molecules adhering 
either through physisorption or the formation of a self- 
assembled monolayer. The nanowires are then broken in 
a controlled manner, either through a mechanical process 
or via electromigration. Ideally, after the breaking pro- 
cess a molecule will be situated in the newly formed gap. 
The broken ends of the wire can then act as source and 
drain electrodes while a capacitively coupled conducting 
substrate can function as a gate electrode (Fig. QJ). 

The tunneling probability decreases exponentially with 
distance between the source and drain, ensuring that in- 
terelectrode conduction is dominated by an extremely 
small volume, comparable to that of a single molecule, 
situated at the point of closest interelectrode separation. 
One can assume that only a fraction of the devices that 
are broken will have an appropriately sized gap contain- 
ing an individual active element. It follows that the for- 
mation of a working sample using any breakjunction tech- 
nique is highly statistical in nature, and many nanowires 
are often required to achieve an operational device. 

In the mechanical breakjunction technique a post is 
positioned against the back of the substrate beneath the 
nanowires, while the substrate edges are rigidly held 
in place. The post is slowly extended by means of a 
piezoelectric stack. The substrate bends as the piezo 
stack is extended, which results in the nanowire being 
stretched until a small break forms (Fig. QJ)).— £ The 
large geometric reduction ratio between the movement 
of the piezo and the nanowire assures a high degree 
of control over the resulting gap size. An alternate 




FIG. 1: a) Schematic of STM measurement setup. b) 
Schematic of mechanically controlled break junction appara- 
tus. Top: SEM image of Au nanowire suspended 40 nm above 
the substrate. (Adapted from Ref. 6 with permission from 
the American Physical Society), c) Idealized single molecule 
break junction device configuration with capacitively coupled 
gate electrode, d) SEM image of lithographically defined Au 
nanowire defined by e-beam lithography on a Si/SiC>2 sub- 
strate. Inset: higher resolution image of similar device after 
electromigration has created a nanometer-size gap. 

method is referred to as the electromigration technique,— 
and has become the most widely utilized means of 
fabricating three-terminal single molecule transistor 
(SMT) devices^— A current is ramped between the 
source and drain electrodes of a nanowire. When the 
power dissipated reaches a certain point a nanometer 
size gap opens up at the narrowest portion of the 
nanowire through the process of electromigration and 
Joule heating (Fig. [TJi). 

BASIC TRANSPORT IN SINGLE MOLECULE 
DEVICES 

Establishing that transport is through a single 
molecule of interest has presented a challenge to ex- 
perimentalists investigating these devices. A STM pro- 
duces a map of the topography and the DOS of a surface 
whereby one can visually search the substrate surface for 
an isolated molecule. A high quality STM may reveal a 
detailed picture of molecular structure sufficient for un- 
equivocally verifying the relevant molecular speciesJ^ 

The situation in break junctions is more difficult. 
Whether or not conduction in a breakjunction device 
is occurring through a single molecule of the desired 
species is usually ascertained through the analysis of 
transport properties, as discussed at length in several 
references.— ~— A number of characteristics indicative of, 
or antithetic to, conduction through a single molecule are 
evaluated before a more in-depth analysis is attempted. 
In addition to determining that charge is passing through 
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only one molecule, one may search for attributes in 
the transport data that uniquely identify the particular 
molecule of interest.— 

The primary tool for characterizing an individual de- 
vice is the measurement of the differential conductance, 
dl/dV, as a function of source-drain bias, V sc i, calcu- 
lated from I — V traces or measured directly using stan- 
dard low frequency lock-in techniques. The dependence 
of this conduction data on temperature, T, and, in three- 
terminal devices, on gate voltage, V g , provides a wealth 
of information about conduction mechanisms. 

It is important to consider the relationship between 
conductance in these solid-state situations and the usual 
discussion of electron transfer theory in chemistry^ - — 
Here the source and drain electrodes are analogous 
to donors and acceptors. Furthermore, in the usual 
solution-based chemistry case, fluctuations in the local 
environment (e.g. rearrangements of polar molecules or 
nearby ions) are critical to the reaction energy landscape, 
whereas in the solid state limit there are no such fluc- 
tuations. Instead, electrostatic interactions between the 
molecule and a solution environment Angstroms away are 
replaced here by interactions with the electrodes. 

The tunneling barriers between the electrodes and the 
active element arc defined by the geometry and chem- 
istry of the molecule-electrode interfaces. The end groups 
used to form chemical bonds to metallic electrodes have 
been a subject of particular interest as they have a di- 
rect impact on the molecule/source and molecule/drain 
coupling strengths, respectively labeled T$ and Tp. In 
an experiment using a modified STM configuration, cur- 
rent was recorded at a fixed bias voltage while repeatedly 
forming and breaking Au point contacts in the presence 
of molecules. Comparing conductance histograms com- 
prised of many conductance traces measured in the pres- 
ence of thiol, isonitrile, and amine terminated aromatics, 
Venkataraman et al. found that the choice of end-group 
greatly affected the distribution and the average values 
of measured off-resonant conductance through the junc- 
tions (Fig. [2^i) J£ More reproducible conductance values 
point to a narrower distribution of end- group/electrode 
bonding geometries. For example, amines are found to 
bind preferentially to undercoordinated Au structures. 

The most elementary models treat the two contacts be- 
tween the chemical end groups and the metal as isolated 
tunnel barriers, with low-energy accessible states in the 
molecule through which sequential tunneling can occur. 
While this can be a useful representation at times, as dis- 
cussed below, such an interpretation would not take into 
account the detailed molecular structure of a particular 
active element, which can play an important role in the 
resulting transport characteristics. For instance, in stud- 
ies of off-resonant conduction, a molecule's lengt h 10 i 19 
(Fig. [2}d) as well as its conformation^ (e.g. the relative 
twist angle between two rings of a 7r-conjugated biphenyl 
molecule) (Fig. 0:,d) exhibit distinct relationships with 
measurements of the average junction conductance. It is 
also generally accepted that in such off-resonant trans- 




FIG. 2: a) Conductance histograms measured as junc- 
tions are pulled apart with 25 mV bias in the presence 
of 1,4-benzenediamine (blue), 1,4-benzenedithiol (red), 1,4- 
benzenediisonitrile (green), and without molecules (yellow) 
shown on a log- log plot. Inset: same data on a linear plot 
with a Gaussian fit to the peak (black). Amine terminated 
molecules result in the most well defined conductance, b) Off- 
resonant conductance is measured through alkanediamines 
with different numbers of methylene groups in the chain. The 
x's denote the center peak position in the conductance his- 
tograms plotted against the number of methylene groups in 
the alkane chain. Junction conductance decreases with an in- 
creasing number of methylene groups. (Reproduced from Ref. 
19 with permission from the American Chemical Society), c) 
Subset of the biphenyl series studied, shown in order of in- 
creasing twist angle or decreasing conjugation. Conductance 
histograms are obtained from measurements through the se- 
ries of molecules with internal twist angles increasing from 
molecule 2 through 8. d) Data points indicate histogram 
peak position for each molecule plotted against cos 2 #, where 
6 is the twist angle calculated for each molecule. The planar 
conformation has the highest conductance. (Reproduced from 
Ref. 20 with permission from the Nature Publishing Group). 



port the derealization of the molecular orbitals make 
conjugated molecules more conducting than typically in- 
sulating saturated molecules. 

In the case of resonant conduction, it is often appro- 
priate to model the system as a double barrier tunneling 
junction, in which the molecule is treated as a quan- 
tum dot and the details of its structure are largely ig- 
nored. In this limit it follows that the electrostatic charg- 
ing energy of a single-molecule device is non-negligible 
due to the small size and correspondingly small capac- 
itance of the active element. An additional electron 
is blocked from tunneling onto the molecule until this 
charging energy, or on-site repulsion ([/), is overcome. 
This suppressed conduction is known as the Coulomb 
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FIG. 3: a) Energy level diagram of double barrier tunneling 
junction, b) Portion of an idealized stability diagram of a 
single molecule device displaying dl/dVsd (shown as relative 
brightness) as a function of V a d and V g . c) Features visible 
in the stability diagram corresponding to: (1) Resonant tun- 
neling, (2) Resonant tunneling through on excited state, (3) 
Elastic cotunneling, (4) Inelastic cotunneling. d) Data from 
a real C6o-based single molecule transistor device exhibiting 
some of these features. Dashed lines indicate blockaded re- 
gions (i.e. Coulomb diamonds). White arrow indicates ex- 
citation corresponding to the 33 meV intra-cage vibrational 
mode of Ceo- 



blockade. The "particle-in-a-box" spacing between dis- 
crete single-particle energy levels in a molecule, A, is also 
non-negligible. At equilibrium the spacing between single 
particle energy states is the energy spacing between the 
highest occupied molecular orbital level (HOMO level) 
and the lowest unoccupied molecule orbital level (LUMO 
level), commonly referred to as the HOMO-LUMO gap 
(Fig. Ki). 

When no molecular state is available between the re- 
spective Fermi levels of the source and drain (fix and /i2), 
there is no resonant tunneling, and current is said to be 
blockaded. At sufficiently large V s d = (hi — ^2)/e, an 
accessible energy level is shifted into the energy range 
between \x\ and /Z2, overcoming the charging energy and 
the HOMO-LUMO gap, and allowing sequential resonant 
tunneling through the molecule. This leads to a peak in 
dl/dV. 

In the limit of comparatively weak molecule- metal elec- 
tronic coupling (Ts, <C A), a capacitivcly coupled 
gate electrode can shift the discrete molecular levels ener- 
getically relative to [i\ and fi2- Thus gate bias, V g , shifts 
the range of V s d over which current is blockaded. Traces 
of dl/dV measured as a function of V s d over a range 
of incrementally changing values of V g may be compiled 
to form a "stability diagram", illustrating where in the 
V s d — V g parameter space, and with what relative mag- 
nitude, conduction is allowed (Fig. [3]b,d). The result- 
ing diamond pattern (Coulomb diamonds) shows regions 



of suppressed tunneling, and therefore regions of stable 
charge states, with fixed average molecular charge. When 
V g brings a discrete level into alignment with fxx = 
resonant tunneling may occur at zero bias. This point 
in the stability diagram where the blockade vanishes is 
referred to as a "charge degeneracy point" . In chem- 
istry language, charge degeneracy points are the V g values 
where redox transitions between molecular charge states 
take place. 

Further analysis of the differential conductance mea- 
surements can be used to uniquely identify the (species 
of) active element. These features include the num- 
ber of charge states (i.e. Coulomb diamonds) accessi- 
ble by application of a gate voltage; the spacing between 
charge state transitions (i.e. distance between adjacent 
Coulomb diamonds) , calculated from the maximum mag- 
nitude of gate bias required to add or remove one elec- 
tron from the molecule; the size of the electron addi- 
tion energy, determined from the maximum magnitude 
of source-drain bias required to overcome the blockaded 
region; the number of characteristic tunneling threshold 
slopes bordering blockaded regions; and the number and 
spacing of energy levels corresponding to excited states. 

Sequential tunneling processes are also possible in 
which the electron tunnels elastically via an excited state 
of the molecule, and in which the electrons loses energy 
due to tunneling events involving excited states of the 
molecule. The signature of this is the presence of addi- 
tional dl/dV peaks that run parallel to the edges of the 
blockaded region in stability diagrams. Excited states 
of a molecule may originate from both local vibrational 
modes and excited electronic states of the system. The 
energy of these excitations can be determined from the 
source-drain voltage at which they intercept the conduc- 
tance gap (Fig. [3}D-d). If these quantized excitations 
can be correlated with known molecular properties, they 
may serve as a fingerprint for molecular identification. 

Higher order tunneling events are another common 
means of transport in single molecule devices. Cotun- 
neling, for instance, may occur in single molecule devices 
at bias voltages that do not overcome the Coulomb re- 
pulsion. This becomes more apparent when the coupling 
between the dot and leads is increased. During elastic 
cotunneling an electron tunnels into and out of the same 
energy level, leaving the molecule in its ground state 
(Fig. [3Jd,c). This is the process that Marcus/Hush the- 
ory refers to as superexchange. In this off-resonant limit 
the molecule acts as an effective tunneling barrier i2I 

Inelastic cotunneling refers to the process by which an 
electron tunnels onto and off of the molecule through 
two different energy levels, leaving the molecule in an 
excited state. Inelastic cotunneling may occur when the 
source-drain bias is such that eV s d — A, where A is the 
energy level spacing between the ground state and the 
first excited state (Fig. [3]b,c)j 22 ' 23 Inelastic cotunneling 
is the process responsible for inelastic electron tunneling 
spectroscopy (IETS)^^ 

In the next section, we discuss another higher order 
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tunneling process involving spin that produces a zero 
bias resonance observed in some single-molecule devices. 
This resonance is associated with the formation of a 
highly correlated electron state known as the Kondo 
effect. 

THE SPIN-i KONDO EFFECT 

At room temperature the electrical resistance of a typ- 
ical metal is dominated by electron-phonon scattering. 
As T is reduced, this resistance decreases since electrons 
can travel more easily as the lattice vibrations of the 
metallic crystal decrease. Neglecting superconductivity, 
as T — > the resistance for most metals is expected 
to saturate to a constant value, with electron scatter- 
ing set by residual impurities and structural defects. It 
was discovered in the 1930s that some metals instead ex- 
hibit an increase in resistance at low temperatures (Fig. 
[4k)3i It was subsequently shown that this low temper- 
ature resistance increase is correlated with the presence 
of magnetic impurities^ In 1964 Jun Kondo presented 
a theory that was able to successfully account for this 
resistance anomaly, which today we know as the Kondo 
effect^ The increased low temperature resistance was 
shown to originate from an increase in the effective scat- 
tering cross-section of the magnetic impurities. A cloud 
of conduction electrons around the impurities serves to 
screen the impurity magnetic moments and increase the 
scattering of electrons near the Fermi level (Fig. [4}d,c). 

The ground state of the (impurity + conduction elec- 
trons) system is a many-body spin singlet. At its root 
the Kondo effect is a consequence of the Pauli exclusion 
principle and strong electron-electron repulsion. For a 
spin- 1/2 impurity, absent interactions, a conduction elec- 
tron could sit on the impurity, pairing up with the spin of 
the local moment, thereby forming a singlet. Including 
interactions, the strong Coulomb repulsion makes this 
energetically prohibitive. However, at low temperature 
higher order tunneling processes may allow the delocal- 
ized conduction electrons to develop a polarization that 
compensates for the localized magnetic moment. The 
quenching of the magnetic moment by these higher order 
processes is what serves to "screen" the unpaired impu- 
rity spin and enhance the scattering of the conduction 
electrons while forming a new many-body state with to- 
tal spin 0. 

The Kondo effect is inherently a many-body phe- 
nomenon that can appear when there is some localized, 
non-zero spin degree of freedom coupled to a sea of 
conduction electrons. Single-particle approaches alone 
cannot account for the existence of this effect. Using 
perturbation methods, Kondo showed that the antifer- 
romagnetic interaction between the Fermi sea and the 
magnetic impurities leads to the observed logarithmic 
rise in electron-impurity scattering as temperature is re- 
duced. While self-consistent density functional theory 
methods have shown accurate predictions for some prop- 
erties of molecular junctions, these ab-initio approaches 




FIG. 4: a) Conductance vs. temperature through a double 
barrier tunneling junction with and without the Kondo ef- 
fect, b) Magnetic impurity in a sea of conduction electrons 
above the Kondo temperature Tk- c) Formation of the Kondo 
screening cloud when T < Tk- Kondo transport mechanism: 
d) initial state, e) virtual state (during which time an electron 
can hop off of the impurity), f) final state. The replacement 
electron of the final state can be of either spin orientation, g) 
Density of states as a function of the equilibrium Fermi level. 
(Adapted from Ref. 38 with permission from The Institute of 
Physics.) 



fail to capture the strong electron correlation effects lead- 
ing to the conductance resonance associated with the 
Kondo statej 29 i 30 A complete physical description re- 
quires non-perturbative methods, namely Renormaliza- 
tion Group (RG).— ~— Numerical RG techniques were 
later developed that could accurately calculate trans- 
port properties, such as resistivity, over a broad range 
of temperatures . 35 ' 36 

An idealization of this system is the Anderson single- 
level impurity model, where the magnetic impurity is 
viewed as having one energy level with an unpaired elec- 
tron at energy level £o (Fig. [4}i)^2i All energy levels 
below £o are fully occupied with electrons, and all en- 
ergy levels above eo are unoccupied. Neither the levels 
above nor below eo contribute to the Kondo resonance. 
When £o is below the equilibrium Fermi level, the un- 
paired electron is trapped with spin- 1/2. It is classically 
forbidden to bring the electron out of the impurity with- 
out adding energy into the system. However, the un- 
certainty principle does allow for this to happen via an 
exchange process whereby the electron is allowed to jump 
off the magnetic impurity and onto an electrode on a time 
scale ~ /i/|eo| (Fig. [4^). During this time, another elec- 
tron from the Fermi sea must jump onto the impurity to 
replace the one that has left (Fig. |4]f ) . The spin of this 
new electron may be either spin-up or spin-down. The 
coherent superposition of many such co-tunneling events 
results in the screening of the local spin, thereby pro- 
ducing the Kondo resonance (Fig. HJg). In a bulk metal 
such a resonance scatters electrons with energies near the 
Fermi level. Because the same electrons are responsi- 
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ble for the low-temperature conductivity, this scattering 
leads to an increase in resistance. In the single-molecule 
device configurations relevant to us, all electrons must 
travel through the active element, and transport is dom- 
inated by tunneling through a single magnetic site. The 
enhanced scattering in this case is forward scattering that 
couples states from the source and drain electrodes, lead- 
ing to a corresponding increase in conductance as T — > 
(Fig. UtOM^i 

The state giving rise to the Kondo effect, consist- 
ing of a localized degenerate state coupled to an elec- 
tron reservoir, is well characterized by the Anderson 
Hamiltonian,— 



H 



Anderson 



ka 



{ v kd ] a c ka + vl<* ka d a ). (1) 




FIG. 5: Two devices used to observe the Kondo effect, a) 
Scanning electron micrograph of an electrostatically defined 
GaAs/AlGaAs quantum dot. b) Atomic force microscope im- 
age of a device consisting of a 2 nm thick nanotube bundle 
with 30 nm thick gold contacts separated by 300 nm. Sketched 
beneath the image are the 300 nm thick Si02 layer and a 
highly p-doped Si substrate, ((a) and (b) reproduced from 
Refs. 44 and 54, respectively, with permission from the Na- 
ture Publishing Company). 



The first two terms account for the electrons in the leads 
and dot, respectively, where ct (cka) creates (destroys) 
an electron in the leads with momentum k, spin a, and 
energy e k ■ The third term represents the Coulomb inter- 
action of the electrons on the island, and the last term 
describes the tunneling between the dot and the leads 
with amplitude v k . 

With a change of variables, the Schrieffer- Wolff unitary 
transformation can be performed on Eq. (JTJ) . Retaining 
terms up to second order in Uk, while eliminating v k to 
first order, one can arrive at the Kondo Hamiltonian, in 
which only the singly-occupied state is allowed on the 
dot. 
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The second term of Eq. ^ describes the interaction 
between the conduction electrons and the local moment 
as an antiferromagnetic coupling ) 28 ' 41 ! 42 



H,. 



J &cond ' S 



molecule ■ 



(3) 



where S is the net spin of the quantum dot, s is sum of 
spin operators of the conduction electrons, and J repre- 
sents the strength of the effective spin-spin interaction. 

The binding energy of the singlet state (with total spin 
= 0) formed between the dclocalized conduction electrons 
and the localized magnetic moment is not simply pro- 
portional to J. The characteristic energy scale associ- 
ated with the (continuous) transition into and out of the 
Kondo state is given by the Haldane relation^ 



-weo(-s + U)/rU 



(4) 



which gives the Kondo temperature for the Anderson 
model, where U is the on-site repulsion and £o is the 
energy level of the molecule through which transport oc- 
curs, as depicted in Fig. [4}1. Thus Tk is exponentially 
dependent on T, the sum of the widths of the tunneling 



barriers, which is in turn extremely sensitive to the pre- 
cise molecule-electrode coupling. We then find that the 
Kondo energy scale (and the associated zero bias reso- 
nance) will be most evident in devices that possess rel- 
atively strong molecule-electrode couplings, which also 
leads to increased overall conductance. 

Renewed interest in Kondo physics took place in the 
late 1990s after the first experimental demonstration of 
the single impurity Kondo effect was accomplished using 
an electrostatically defined semiconductor quantum dot 
(QD), in which the conduction electrons interact with a 
localized "puddle" of electrons comprising the magnetic 
impurity (Fig. [5fc.) i 44 ' 45 A zero bias Kondo resonance 
was observed due to the screening of the local moment 
by the conduction electrons. At around the same time, 
in 1998 evidence of the Kondo effect was seen in conduc- 
tance measurements through single magnetic adatoms on 
a conducting surface using a scanning tunneling micro- 
scope (STM) 

The general behavior of the zero bias conductance am- 
plitude in a quantum dot can be shown to obey a log- 
arithmic temperature dependence. Goldhabcr-Gordon 
et al. developed an empirical form derived from a fit 
to the renormalization group analysis ; 32 ' 46 which more 
accurately captures the temperature dependence of the 
equilibrium conductance, and can be written as 



G(T,0) 



(i + (2V* - i){T/T K yy 



+ G b 



(5) 



where s = 0.22 for a spin-1/2 impurity, Go is the 
peak conductance minus G& in the limit that T — > OK, 
Go = G(0, 0) — Gfc, and Tk is defined as the value of 
T for which G(T,0) - G b = G /2 (Fig. [BJ. Note that 
in STM and single-molecule junction experiments, there 
can be an additional "background" contribution to the 
conductance, G^, due to non-resonant processes such as 
direct tunneling between the source and drain electrodes. 
The relative source-drain coupling can be deduced from 
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FIG. 6: a) Traces of dl/dV vs. V s d for several temperatures 
ranging from 1.8 K (tallest/black) to 53 K (shortest/maroon) 
for a single molecule device constructed via the electromi- 
gration technique. Amplitude of Kondo resonance decreases 
with increasing temperature, b) Fit of equilibrium conduc- 
tance vs. temperature, normalized with respect to Go and 
Tk, respectively, using Eq. \5\ yielding Tk ~ 52 K. (Adapted 
from Ref. 50). 




FIG. 7: a) Magnetic-field dependence of the Kondo reso- 
nance. (Reproduced from Ref. 9 with permission from the 
Nature Publishing Company), b) Peak positions as a func- 
tion of applied field (blue x's). The peak splitting varies 
linearly with magnetic field greater than a critical value . 52 ' 53 
Dashed lines represents the predicted Zeeman splitting. 



the magnitude of the resonance peak as the temperature 
approaches K using the equation 4 ^ 



Go — Gb 



2e 2 (4r s r D ) 
h (r s + r D ) : 



■f(T/T K ), 



(6) 



where f(T/T K ) = (1 + (2 1 /* - 1)(T/T K ) 2 )~ S . 

Strong correlations between conduction electrons and 
a single Ti atom on a Ag(100) surface leading to the 
Kondo effect were observed using a STM apparatus and 
used to study the temperature dependent broadening of 
the associated zero bias resonance.— Approximating the 
Kondo resonance peak as a Lorentzian, the full width at 
half max (FWHM) has a temperature dependence that 
can be approximated as 



FWHM = - y /( 7T k B Ty + 2{k B T K f, 
e 



(7) 



where ks is the Boltzmann constant. The system pa- 
rameters, Tk and Go, can then be extracted using Eqs. 
(O and ([7]). As the temperature is increased, higher or- 
der transport processes become relevant, and deviations 
from these forms are expected. The enhanced V s d = 
V conductance of the Kondo resonance rapidly decreases 
as a function of applied bias. Quantum interference be- 
tween Kondo processes and non-Kondo processes {e.g. 
direct tunneling) can modify the observed resonance in 
dl/dV vs. V s d, resulting in a Fano lineshape. Depending 
on the relative complex amplitudes for the processes in- 
volved, these lineshapes can take the form of peaks, dips, 
or intermediate asymmetric structures ! 48 ' 49 

The equilibrium conductance of the resonance peak 
can be suppressed by an applied magnetic field which lifts 
the degeneracy of the up and down spin states through 
which transport occurs. In the presence of an applied B- 
field, the enhanced conductance can be recovered through 
the application of a finite bias thereby revealing a split- 
ting of the original Kondo resonance peak (Fig. [7J . The 
split peaks are expected to separate by an amount equal 
to twice the Zeeman energ y 10 i 51 for applied fields larger 



than a critical value, B c ~ 0.5Tk~2. These descriptions 
of the temperature and S-field dependence of the Kondo 
resonance are commonly used to verify a resonance peak 
originating from a single spin-1/2 impurity. 

The energy scales relevant to the Kondo effect in a 
semiconductor QD (i.e. Tk and A) can be prohibitively 
low, making observation of this phenomenon difficult to 
achieve. Common values of Tk in these structures, for 
instance are in the few hundred milliKelvin range. This 
is caused by the relatively weak confinement of the elec- 
trons in the dot, which produces a larger effective dot 
(i.e. impurity) size compared to what is found in a sin- 
gle molecule device. The larger size of the dot leads to 
a correspondingly smaller charging energy, U, which by 
the Haldane relation, Eq. Q, results in a smaller Tk- 
The single impurity Kondo effect can also bc observed in 
carbon nanotubes devices (Fig. (5]b) ) 54 ' 55 and SMTs£~— 
when an unpaired electron exists on the nanotube or 
molecule, respectively, at temperatures below the mean 
energy level spacing of the device. In carbon nanotubes 
the relatively greater confinement, compared to that of a 
semiconductor QD, has led to observed Kondo temper- 
atures typically around a few Kelvin. Single molecule- 
based devices most often exhibit Kondo temperatures 
in the tens to hundreds of Kelvin range. Tk in single 
molecule devices exceeding 100 K was first reported by 
Yu et at in 2004 using Ceo molecules^ 1 - This large mag- 
nitude of Tk is more readily accessible, making molecules 
a desirable system for studying Kondo physics. 

Kondo resonances in molecular devices were first ob- 
served in some of the earliest studies of transport in single 
molecule break junction structures. Park et at demon- 
strated conduction through two different molecules 
([Co(tpy-(CH 2 )5-SH) 2 ] 2 + and [Co(tpy-SH) 2 ] 2 +), both 
of which consisted of a Co-containing coordination 
complex between two terpyridinyl linker molecules with 
thiol end groups.— Using these molecules they were 
able to observe examples of Coulomb blockade and the 
single impurity Kondo effect. Liang et al. fabricated 
a single molecule brcakjunction device with individual 



8 



divanadium (V 2 ) molecules ([(N,N',N"-trimethyl-l,4,7- 
triazacyclononane)2-V2(CN)4(/i-C4N 4 )])r£ and were 
able to demonstrate charge state transitions by gate 
modulation, exhibiting the evolution of transport be- 
tween the Coulomb blockade regime and the Kondo 
regime. Since 2002 Kondo resonances have been ob- 
served in molecular devices involving a wide variety of 
molecular active elements, as reported in the literature. 




SELECTIVE MANIPULATION 
KONDO RESONANCE 



OF THE 



Selectively modifying single molecule devices in order 
to affect changes in the observed zero bias resonance can 
serve as a means of probing Kondo physics in these sys- 
tems. Experimentalists have been able to manipulate 
both the molecular active element as well as its local en- 
vironment through their choice of materials and by more 
direct physical probes. The resulting impact to the reso- 
nance peak has allowed for the investigation of the man- 
ner in which Kondo correlations arc influenced by these 
factors in addition to associated phenomena. 

The zero bias resonance peak is expected to split 
by approximately 2g[isB in the presence of an applied 
magnetic field, as stated previously; however, compet- 
ing magnetic exchange interactions can significantly al- 
ter this response, as demonstrated in an experiment by 
Pasupathy et al. Cgo-based single molecule breakjunc- 
tion devices were fabricated using ferromagnetic nickel 
electrodes (Fig. [8^0 The source and drain electrodes 
were designed to have different magnetic anisotropies, in 
order to ensure that they undergo magnetic reversal at 
different values of B. Although the Kondo effect can 
be suppressed due to a competing ferromagnetic ground 
state, zero bias resonances were still achieved in some of 
these devices in which the molecule is strongly coupled to 
the electrodes, but the resonances were often found with 
peak splittings too large to be accounted for by a local 
field produced by the Ni electrodes. 

It was found that the interaction between the spin po- 
larized electrodes and the molecule results in a large local 
exchange field (> 50 T). The spin asymmetry in the cou- 
pling between the molecule and electrodes, due to quan- 
tum charge fluctuations, leads to a spin-dependent renor- 
malization of the molecular energy levels, breaking the 
spin degeneracy, and results in a splitting of the reso- 
nance peak. Furthermore, the splitting of the zero bias 
resonance is strongly dependent upon the relative ori- 
entation of the electrode magnetizations, which can be 
controlled with a small external B-ficld (< 100 mT) (Fig. 
[8)d,c). The manipulation of electrode magnetizations is 
confirmed by measuring the junction magnetoresistance 
(JMR) of molecule-free electromigrated Ni junctions. For 
transport in the Kondo resonance, the JMR can be am- 
plified by turning the electrode magnetization from par- 
allel to antiparallel, since the Kondo resonance occurs 
closer to the Fermi energy when the magnetization of the 
leads is antiparallel. In this respect, the Kondo resonance 
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FIG. 8: a) False color SEM image of break junction fab- 
ricated with Ni electrodes. Inset: High resolution image af- 
ter electromigration. b) and c) Differential conductance vs. 
source-drain bias for two Ceo-based devices with Ni electrodes 
at T = 1.5 K. Blue (lower amplitude) lines: electrode mag- 
netizations parallel, (b) B = -310 mT, (c) B = -250 mT. 
Green (higher amplitude) lines: electrode magnetizations ap- 
proximately antiparallel, (b) B — -10 mT, (c) B — -15 mT. 
(Reproduced from Ref. 56 with permission from the AAAS) 



serves as an indicator when probing the relative orienta- 
tion of the source and drain electrode magnetizations and 
their interaction with the local magnetic moment of the 
molecule. 

The Kondo effect can also be used to probe the in- 
teraction between the unpaired spin of the local mo- 
ment and magnetic impurities in the adjacent leads. 
The Rudcrman-Kittel-Kasuya-Yoshida (RKKY) interac- 
tion describes the spin-spin interaction by which mag- 
netic impurities embedded in an electron sea can inter- 
act with each other j££ Elastic spin-flip processes leading 
to the Kondo resonance may be suppressed by compe- 
tition from the RKKY exchange interaction, which may 
essentially freeze the spin on the impurity, thereby reduc- 
ing the conductance near V sc t = V and leading to an 
observed peak splitting. 

An investigation by Heersche et al. observed Kondo 
resonances in gold break junction devices resulting from 
an unpaired spin on a small gold grain formed during the 
electromigration procedure^ A small number of cobalt 
impurities were intentionally deposited in the electrodes 
during the evaporation step of the fabrication process. 
They concluded that the (RKKY) interaction between 
the fixed magnetic impurities in the electrodes and an 
unpaired spin on the gold grain, which serves as the active 
clement, cause the zero bias Kondo resonance to split. 

The magnetic-field dependence of the spacing between 
the split peaks depends on the sign of the associated 
RKKY interaction, /. The interaction can be either fer- 
romagnetic (/ < 0) or antifcrromagnctic (/ > 0), however 
both cases will generally suppress the s — 1/2 Kondo ef- 
fect. Depending on the type of interaction, the net spin of 
the active element and the impurity spins together can 
form a triplet (S = 1) state or an underscreened spin- 
1 state. Because the energy between the singlet ground 
state and the triplet state decreases with \B\, an external 
field can restore the Kondo effect if the RKKY interaction 
is antiferromagnetic j ' The zero bias resonance peak is 
restored at B = I/gfj,sB, where the singlet and triplet 
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FIG. 9: Magnetic-field dependence of the Kondo resonance 
in a nanowire break junction with magnetic impurities em- 
bedded in the leads, a) The singlet-triplet transition energy 
(vertical arrow) decreases, then increases with B-field for an 
antiferromagnetic interaction. The triplet-singlet energy al- 
ways increases with field for a ferromagnetic interaction, b), 
c) Traces of conductance vs. applied bias obtained at differ- 
ent values of the external magnetic field, increasing from B 
= T (bottom line) to B = 9 T (top line) in steps of 0.9 T. 
Data taken at T = 250 mK (b) Resonance peak restored at fi- 
nite field for antiferromagnetic interaction between two spins. 
Traces offset for clarity, (c) The peak separation increases 
linearly with \B\ for a ferromagnetic interaction. Inset: The 
peak separation at B = T is determined by extrapolating 
from peak positions to zero field and yields 1.4 ± 0.3 meV. 
(Reproduced from Ref. 57 with permission from the Ameri- 
can Physical Society). 



states are degenerate and the external field compensates 
the antiferromagnetic interaction. In the case of a fer- 
romagnetic interaction the peak splitting will increase 
monotonically with \B\ because the energy between the 
triplet and singlet states also increases (Fig. [9j. 

The magnetic field dependence can therefore be used 
as a probe to distinguish between ferromagnetic and an- 
tiferromagnetic interactions between the local moment 
and conduction electrons. While the authors have used 
an entirely metallic system, the local moment could be 
equally well represented in future experiments by an un- 
paired spin on a molecule. The break junction system 
used in this experiment represents a flexible platform to 
study the interaction between the spin of the local mo- 
ment and adjacent static magnetic impurities. 

Controlled physical alterations in the structure of sin- 
gle molecule devices provides another means of probing 
Kondo correlations. This has been accomplished by both 
mechanical and electrical means capable of changing the 
precise arrangement of the source and drain electrodes, 



as well as the detailed structure and conformation of the 
active element, and its interaction with the local envi- 
ronment. Studies of this nature have been performed 
predominantly with the STM configuration, however the 
use of a mechanically controlled break junction has also 
been employed. 

Single molecule devices constructed with a STM may 
not present a scalable device concept; however, they do 
provide a unique platform for probing the interaction of 
conduction electrons with molecular orbitals and individ- 
ual magnetic ions. A STM possesses the unique ability 
to both image the DOS of individual molecules on a con- 
ducting surface, and to manipulate the detailed chem- 
ical environment via mechanical and electrical means. 
When the metallic tip is positioned over a magnetic ion 
or molecule, Kondo exchange interactions may result in 
a resonance at the equilibrium Fermi level observed in 
measurements of dl/dV. The characteristic energy scale, 
Tk, extracted from the measured Kondo resonance indi- 
cates the relative strength of the coupling between the 
local magnetic moment and the delocalized conduction 
electrons. Repeating these measurements at many points 
will provide a spatial map of the Kondo resonance, there- 
fore revealing an image of this coupling strength over the 
surface of molecules. 

Several studies of Kondo correlations have exploited 
the ability of an STM to exert control of the precise sur- 
face chemistry and conformation of a single molecule or 
small group of molecules. These studies have focused on 
Co-containing molecules adsorbed on nonmagnetic con- 
ducting substrates. The authors of Zhao et al. used 
a STM to measure differential conductance through a 
cobalt phthalocyaninc molecule adsorbed on a Au(lll) 
surface-^ Initial measurements exhibited no zero bias 
resonance since the conduction electrons did not inter- 
act strongly with the Co ion. Dehydrogenation of the 
ligands {i.e. cutting off the hydrogen atoms with volt- 
age pulses from the STM tip) allows the four outward 
extending orbitals of this molecule to chemically bond 
to the substrate. This bonding changes the conforma- 
tion of the molecule and its coupling to the substrate, 
and hence its interaction with the conduction electrons. 
After this conformational change, a clear Kondo reso- 
nance can be observed near the Fermi surface. Similar 
work was performed by Wahl et al. in which measure- 
ments of Kondo coupling in cobalt tetracarbonyl were 
analyzed before and after tip-induced dissociation of the 
molecule.— 

A related experiment investigated the TBrPP-Co 
(5,10,15 , 20-tetrakis (4-bromophcnyl) porphyrin-Co) 
molecule, which is found in two conformations, "planar" 
and "saddle", both of which anchor themselves onto 
a Cu(lll) surface via their four bromine atoms (Fig. 
110b . b, and c). A single voltage pulse of 2.2 V from 
a STM tip positioned above one of these molecules 
has the ability to switch the complex between the two 
conformations.— A Kondo resonance can be observed 
in molecular devices made with both conformations, 
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FIG. 10: a) STM images of saddle conformation (width 11 
A, length 18 A) and b) planar conformation (length 15.5 A) 
of TBrPP-Co on Cu(lll). c) Corresponding models of sad- 
dle (left) and planar (right) conformation. Blue and pink 
balls represent bromine and carbon atoms, respectively. The 
black ovals indicate the areas providing higher current in the 
STM images. The length of the molecule along the diagonal 
(dashed green arrow) is unchanged upon switching between 
the two conformations. Traces of dl /dV vs. applied bias ex- 
hibiting Kondo signatures for d) a saddle TBrPP-Co molecule 
and e) a planar TBrPP-Co molecule. The solid lines repre- 
sents the Fano line-shape fits to the data. The spectra for 
the saddle and planar conformations have been vertically off- 
set by 0.07, 0.04, respectively. Measurements were taken at 
the center of each molecule. (Reproduced from Ref. 63 with 
permission from the American Chemical Society) 



but they possess distinctively different values of Tk, 
and hence experience different interactions with the 
substrate. The Kondo temperature is decreased in the 
case of the saddle conformation as compared to the 
planar TBrPPCo (Fig. [TO}!, e). 

The authors suggest that the measured Kondo reso- 
nances in these molecules differ because of the potential 
for conduction to originate from more than one path. 
They indicate that a resonance may be generated by 
spin-electron coupling to the substrate through bonding 
orbitals (via the bromine atoms) and by direct coupling 
to the Co ion. When in the saddle conformation, the 
Co atom is lifted away from the surface as the porphyrin 
plane is bent upward, reducing the coupling between the 
Co atom and the free electrons of the Cu(lll) substrate 
(Fig. 110b ). Coupling via molecular orbital bonding is 
therefore likely to be the only conduction pathway avail- 
able in this state. For the planar (higher Tk) conforma- 
tion the Kondo effect may result from the orbital bonding 
as well as Co-surface interaction. This single molecule 
switching mechanism reveals a means of altering the rel- 
evant interactions without destroying or changing the 
chemical composition. 

A similar study indicated that the chemical envi- 
ronment near a molecule also plays an important role 
with respect to Kondo correlations. Beginning with a 
2D molecular assembly of TBrPP-Co, where the cen- 
ter molecule is initially surrounded by six adjacent 
molecules^ Iancu et al. use a STM tip to measure 



the Kondo resonance, and thus Tk, above the center 
molecule. One nearest neighbor (nn) molecule at a time is 
removed from the assembly, after which Tk is again mea- 
sured above the molecule originally at the center (Fig. 
1111) . Measurements of Tk before decomposition of the 
assembly indicate values that are the same as those mea- 
sured inside a full ribbon structure. 

As the number of nn molecules is reduced, Tk is in- 
creased in a controlled manner from 105 K to 170 K. It is 
found that manipulating the nearest neighbor configura- 
tion has a direct effect on the coupling strength between 
the local moment and conduction electrons responsible 
for the Kondo resonance. The molecule located inside the 
cluster will have a reduced interaction with the delocal- 
ized surface state electrons as compared to the molecules 
around the edges. This again represents a means of af- 
fecting the interactions responsible for the Kondo effect 
without affecting the chemical composition. 

A mechanically controlled break-junction can also 
selectively manipulate the molecule-electrode coupling. 
The controlled movement of a piezo actuator is used to 
stretch a nanowire until a small gap appears. In the 
case of J. J. Parks et al& this technique was employed to 
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FIG. 11: a) Sequence of STM images (left), created by remov- 
ing one nn molecule at a time with the STM tip, and corre- 
sponding models (middle). The calculated electron standing 
wave patterns (right) reveal an increasing exposure of the cen- 
ter molecule (red circles) to the surface state electrons. White 
and black colors in the calculated images represent higher and 
lower electron densities, respectively, b) The dl /dV spectra 
are measured at each step by positioning the tip above the 
center molecule (red dots in (a)). The spectra are offset for 
clarity. Horizontal displacements of 3 to 10 meV are taken for 
the spectra representing nn6, nn5, nn4, nn2, and nnl. (c) Tk 
plotted as a function of the number of nearest neighbors. (Re- 
produced from Ref. 64 with permission from the American 
Physical Society) 
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FIG. 12: a) dl/dV traces at various electrode spacings and 
T = 1.6 K for a single mechanically controlled break junction 
device. Smallest spacing corresponds to tallest peak. Largest 
spacing corresponds to shortest peak, b) The Kondo temper- 
ature, Tk, and the relative coupling, Tm/Tl, as a function 
of electrode spacing for same device at T = 1.6 K. The un- 
certainty in the determination of Tk values is < 0.4 K. c) 
dl/dV for another device at T = 1.6 K, exhibiting satellite 
peaks near ±33 mV. Left inset: Schematic of intracage vibra- 
tional mode, d) d 2 I/dV 2 as a function of bias voltage and 
electrode spacing for device in (c). (Reproduced from Ref. 6 
with permission from the American Physical Society). 

produce a Cgo-based single molecule device exhibiting a 
spin-1/2 Kondo resonance, after which the interelectrode 
gap size was further modified (with picometer precision) 
to change the geometry of the molecule-electrode assem- 
bly. Initial values of Tk are extracted from measurement 
of the equilibrium conductance as a function of temper- 
ature. Increasing the spacing between source and drain 
electrodes is shown to tune both the height and width 
of the Kondo resonance, effectively decreasing Tk by as 
much as 37% in one example. The mechanical manipu- 
lation thus has a direct impact on the Kondo exchange 
correlations by altering the relative coupling strength be- 
tween the molecule and the two electrodes. A standard 
theoretical model is able to characterize the changes in 
molecular orbital-electrode interaction (and thus Tk) due 
to electrode displacement. 

Excited states, manifested as finite-bias peaks in 
dl/dV, appeared in some devices at energies consistent 
with the 33-35 mV intracage vibration for Cgo, as in 
Refs. 8 and 11. Similar vibrational sidebands have 
recently been observed in STM measurements of Kondo 
resonances in organic charge-transfer assemblies^ The 
mechanical motion also affects the energy of these active 
vibrational modes, but not in accord with predictions 
based on a simple semiempirical Hamiltonian. This 
suggests that basic approaches used to capture transport 
properties in these devices may be insufficient, and that 
new models or approaches will ultimately be required to 



provide a complete explanation. Further observations 
of transport in the Kondo regime of single molecule 
devices, as discussed in the following section, have led to 
similar conclusions. 

OBSERVATIONS AND COMPARISONS TO 
STANDARD MODELS 

The basic theoretical descriptions of the Kondo state, 
outlined above, have enjoyed a good deal of success in 
their agreement with measured transport properties in 
experiments using semiconductor QDs. This, together 
with the advantage of tunable tunneling barriers and a 
high level of reproducibility, has made electrostatically 
defined semiconductor QDs a canonical system for study- 
ing the Kondo state. Transport measurements in QDs 
made in materials like GaAs have provided a strong foun- 
dation for researchers investigating analogous effects in 
molecular devices. Experimental investigations of Kondo 
correlations in devices where a single molecule is utilized 
as a type of "natural" QD reveal transport properties 
for which the models described earlier cannot fully ac- 
count. These results indicate that sophisticated analysis 
techniques together with a more thorough understanding 
of the underlying physics are required to resolve these 
molecular complications. 

One discrepancy between physical descriptions suffi- 
cient for semiconductor QDs and SMTs can be seen in 
the functional form of the low energy nonequilibrium con- 
duction. Scaling can be seen in the response of the Kondo 
system to different perturbations. Applied bias, tempera- 
ture, and magnetic field may perturb the low energy con- 
ductance of a Kondo resonance in qualitatively different 
ways, but the exponents of their lowest order responses 
are identical. Theoretical treatments suggest that the low 
energy region of every Kondo system can be described 
by the same universal scaling function (i.e. conductance 
collapses onto a single curve) when it is normalized by 
the characteristic energy scale, k B TK^- Limited ex- 
perimental work with GaAs QDs verify that the non- 
equilibrium conductance of the single channel spin-1/2 
Kondo resonance at low bias and temperature adheres to 
a scaling function characterized by two parameters^ 

We recently tested the ability of a single scaling func- 
tion to describe transport in the spin-1/2 Kondo regime 
of single molecule devices. Differential conductance was 
analyzed as a function of both temperature and applied 
bias for a large number of samples (29) fabricated on a 
Si/Si02 substrate. The electromigration technique was 
employed to create nanometer-size gaps in Au nanowircs 
while alternately using two different molecular active el- 
ements: Ceo and a Cu-containing version of the transi- 
tion metal complex in Fig. I13h (bis(2,5-di-[2]pyridyl-3,4- 
dithiocyanto-pyrrolate)Cu(II)). The Ceo molecules can 
be strongly bound to the Au electrodes by charge trans- 
fer from the metal surface to the molecule^ The con- 
jugated complexes self-assemble on the Au nanowires in 
tctrahydrofuran (THF) through loss of the _ CN moieties 
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FIG. 13: a) Assembly sequence on gold and final assembled structure of the Co-containing complex. (From Ref. 14) b) Cyclic 
voltammograms showing reversible change in Co ion charge state, c) Stability diagram of a device containing a single complex 
exhibiting a gate induced transition between the Coulomb blockade regime and the Kondo regime (T = 5 K). White corresponds 
to dl/dV = 3 x 1CP 6 S. The zero bias peak in the right-hand charge state disappears at the charge degeneracy point. (From 
Ref. 71). 



and formation of Au-S covalent bonds (Fig. I13b )i 70 ' 71 
These devices were found to have a variety of Kondo 
temperatures, ranging from 35 K to 155 K, based on the 
temperature-dependence of their zero bias conductance. 
They were also found to possess a range of asymmetries, 
as inferred from the magnitude of the resonance peak as 
T -> K, according to Eq. 

Transport measurements of the temperature and bias 
dependence of the non-equilibrium conductance are in 
agreement with the functional form of the expected scal- 
ing function. Namely, a quadratic power law in tem- 
perature and bias characterized by two scaling coeffi- 
cients. However, values of the extracted scaling coeffi- 
cients are found which are consistent throughout the data 
set, but differ systematically from both theoretically pre- 
dicted values and coefficient values extracted from mea- 
surements of nonequilibrium transport through a GaAs 
dot i 50 i 68 The low energy approximation of the scaling 
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FIG. 14: a) Conductance as a function of applied bias for 
T/Tk < 1.0. A background contribution to the conductance, 
Gb, has been subtracted off. b) Scaled conductance, AG/av 
vs. (eV I 'IcbTk) 2 for four samples, including the data from 
device in (a), using the average extracted values of the scaling 
coefficients, where AG = (1 - G(T, V)/G(T, 0)). The solid 
black line represents the associated universal curve. The color 
bar in (b) pertains to both figures. (Adapted from ref. 50). 



function can be written a o 48 i 50 ' 67 i 68 ' 72 
G(T, 0) — G(T, V) 



c T G 



eV 



T V 



eV 



(8) 



where a and 7 are the scaling coefficients and ct ~ 4.92 
is a constant. We calculated average values of a = 
0.051 ± 0.01 and 7 = 0.107 ± 0.027. This is in contrast to 
values of a = 0.10 ± 0.015 and 7 = 0.5 ± 0.1, for a GaAs 
QD. The consistently smaller values of a and 7 found 
for all of these molecular devices, compared to the GaAs 
case as well as predictions based on the Anderso n 36 ' 73 1 74 
and Kondo 42 ' 67 ' 72 ' 75 models, implies that for a given Tk 
inferred from the temperature dependence of the equi- 
librium conductance, the resonance peaks observed here 
were broader in voltage, and evolved more slowly with 
temperature. 

The uniformity of extracted coefficient values for SMTs 
made with the two different molecules suggests that the 
origins of this discrepancy likely do not depend on the 
detailed molecular structure of the active element. The 
size of the on-site repulsion, U, of the molecule relative 
to the single-particle level spacing, A, may play a role. 
Generic vibrational modes in SMTs, such as a center-of- 
mass oscillation, may also have an effect by modifying the 
temperature dependence leading to a different effective 
Tk than is found in the bare Kondo system, as discussed 
below. 

Further experiments using the same transition metal 
complex as above, alternating between versions of the 
complex containing a Cu and a Co ion, serve to high- 
light some of the distinguishing transport characteristics 
of the single molecule system. Again, single molecule 
breakjunction devices are fabricated on a Si/Si02 sub- 
strate using the electromigration technique. Zero bias 
resonance peaks were observed to have a temperature 
and magnetic field dependence consistent with the spin- 
1/2 single impurity Kondo effect. The start (or end) of 
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the resonance, with respect to a gate bias, was always 
coincident with the charge degeneracy point (Fig. 113b ). 
and similar features were not seen in any control devices. 

Evolution of the Kondo resonance with gate bias was 
found to be inconsistent with the simple model, which 
works well for analogous data obtained in semiconductor 
QDs and carbon nanotube devices. More specifically, it 
was found that Tk{Vq) was distinctly weaker than that 
expected according to Eq. (j4])^ It has been suggested 
that the specific orientation of the active element relative 
to the detailed electrode geometry can result in poor gate 
coupling prohibiting effective gate modulation^! Anoma- 
lously weak gate dependence may also be traced back to a 
number of potentially incorrect assumptions made when 
invoking Eqn. An invalid value of Sq/Y, which is in- 
ferred by normalizing V g by the width of the Coulomb 
blockade charge degeneracy point, may be due to some 
mechanism intrinsic to the SMT system. Screening cor- 
relations in the mixed valence regime can renormalizc the 
measured Y to a value different than the Y relevant to the 
Kondo temperature. 

The energy of intramolecular exchanges may be an- 
other factor requiring further consideration. This energy 
is typically ignored, yet may be comparable to U and A in 
single molecule devices. Electronic structure calculations 
confirm that the coupling to intramolecular excitations 
is strong in both the Co and Cu-containing complexes, 
complicating the resulting molecular orbital degeneracies 
and/or splittings. Ligand field effects, for example, may 
split the degeneracy of the d states, and derealization 
of the majority spin states, extending into the ligands, 
is consistent with large Y values for these systems. De- 
generacies for the unpaired spin in the molecular orbitals 
hybridized with the metallic electrodes may enhance Tk, 
making the normalization of the abscissa correspondingly 
too large.— Intramolecular exchange and the coupling of 
vibrational modes to charge may both increase Tk and 
decrease its gate dependence compared to an identical 
system with no vibrational coupling i 79 ' 80 Among the pos- 
sible influences discussed here, a considerable amount of 
attention has been paid to the impact of vibrational ef- 
fects. 

Local vibrational modes are a distinctive feature in the 
single molecule system distinguishing it from the other- 
wise analogous case of the semiconductor quantum dot. 
In semiconductor dots, optical phonon modes are con- 
siderably higher in energy than the Coulomb charging 
scale, and arc therefore typically ignored in theoretical 
models of the low energy dot properties. In the molec- 
ular case, however, local vibrational modes can have en- 
ergies ~tens of meV, often considerably lower than the 
Coulomb scale. Molecules can move and distort upon the 
addition or removal of an electron, leading to phonon as- 
sisted tunneling, which is commonly observed in stability 
diagrams plotting conductance in the Coulomb blockade 
regime. Conductance features related to intramolecular 
vibration o 6 ' have also shown up in the Kondo regime, 
and similar excitations expose other unresolved transport 



phenomena. 

Inelastic electron tunneling spectroscopy (IETS) re- 
veals features in a portion of the SMTs studied, which 
correspond energetically to vibrational excitations of 
the molecule, as determined by Raman and infrared 
spectroscopy. The tunablity of these systems allows 
for the examination of such features across successive 
charge states. Inelastic cotunncling processes can pro- 
duce peaks, dips, or intermediate structures in d 2 I/dV 2 , 
which are observed in the blockaded regime of an even 
charge state and persist into the Kondo regime in the 
form of satellite features paralleling the zero bias reso- 
nance (Fig. [3)3,0) f 7 ^ These features are a signature of 
IETS, and can be seen to undergo a clear change in line- 
shape and intensity after passing the charge degeneracy 
point (Fig. 115)) . indicating a complicated conduction pro- 
cess and non-trivial transition between the two conduc- 
tion regimes. 

Recent analytical and numerical investigations of 
transport in single molecule devices have considered the 
role of electron-phonon coupling and its relationship with 
Kondo correlations. These effects can be modeled with 
appropriate modifications to the Anderson Hamiltonian. 
A term may be added to Eq. |T]) to account for the local 
vibrational modes of the molecule, 

l 

which arc represented here as a discrete number of har- 
monic oscillator modes. The creation and annihilation 
operators a\ and a l can be used to define a displacement 
operator for a given mode: 

ar, = ./-A- (aj + a,). (10) 

Here mi is an effective mass characterizing the local vi- 
brational mode of frequency W;. Another term, utilizing 
the displacement operator, is employed to describe the 
coupling between these local modes and the molecular 
energy levels, 

Ira 

Here A is a dimensionlcss matrix of coupling constants 
that links molecular level v with local vibrational mode 
I. 

The local displacement of the molecule due to vibra- 
tions, including nanomechanical oscillations, can mod- 
ulate the molecule-electrode tunneling barriers, so the 
portion of the Hamiltonian that describes the tunneling 
of electrons between these components may be altered to 
incorporate an explicit dependence on Eq. (|10[) . 81 Terms 
are also included to account for bulk vibrational modes 
(phonons) of the electrodes 

H ph = Y,m t b\b l , (12) 

9 
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and their coupling to local vibrational modes of the 
molecule, 

H vph =^P lq x l X q . (13) 

lq 

In order to make the resulting Hamiltonian tractable, 
a number of these treatments consider the coupling of 
vibrations to only a single electronic level, and have 
utilized a modified version of the Anderson Hamilto- 
nian (the Anderson-Holstein model) that incorporates 
one vibrational mode of the molecule linearly coupled to 
charge fluctuationsi^Sr— In this approximation, summa- 
tions over the v and I indicies are absent. These studies 
confirm that the nature of the Kondo effect can be mod- 
ified by electron-phonon interactions, and the resulting 
effects will depend on the symmetry and frequency of 
the phonon mode and on the strength of the coupling. 

In the regime of weak electron-phonon coupling (U 3> 
2\ 2 /huio) the same basic Kondo effect is expected to oc- 
cur for a spin- 1/2 impurity, but the interaction poten- 
tially results in an increase of the spin fluctuation energy 
and of the amplitude of charge fluctuations. In the case 
of equal source-drain coupling, the effective charging en- 
ergy is reduced to a new value, 

U eff ^U-2X 2 hoj Q . (14) 

The calculated effect also leads to an increase in the value 
of Tk and a rcnormalization of the electronic energy level, 
Eq, and the single-particle level width, T. These changes 
have been shown to correlate to an altered Tk{V 9 ) and 
a weakening of the influence from an applied gate volt- 
age, as observed in Ref. 76. Another change in gate 
dependence may manifest as an asymmetry in the curve 
of G(V g ) about its maximum^ 

In the strong electron-phonon coupling regime (U <C 
2A 2 /fiwo) the properties of the device are described in 
terms of a charge analog of the Kondo effect i 80 ' 86 In 
this regime the transport characteristics of the molecular 
junction are controlled by the large charge polarizability 
of the molecule. The onsite repulsion, U, can be renor- 
malized to a negative value. This means that the effective 
U produces an attractive interaction between the elec- 
trons, which is described as a polaronic effect and thus 
leads to qualitatively different physics. One consequence, 
for example, is a general insensitivity of the conductance 
to an externally applied magnetic fields More sophisti- 
cated modeling techniques, such as electronic structure 
calculations involving molecules bound to realistic elec- 
trodes while including complex correlations, would also 
greatly facilitate this research. 

Further experimental data will be needed to provide 
support to existing theories and insights into the distinc- 
tive nature of transport in molecular electronic devices. 
For example, an analysis of shot noise in molecular de- 
vices would enable a new method of studying exchange 
correlations, and may serve to contrast fundamental dif- 
ferences in transport behavior between single molecule 




FIG. 15: Colormaps of d 2 I/dV 3 2 d as a function of V a d and 
V g at 5 K for two devices. Smoothing window in V a d is 5 
mV. Brightness scales are -8 x 10" 5 A/V 2 (black) to 3 x 
1(T 5 A/V 2 (white), and -2 x 10~ 5 A/V 2 (black) to 2 x 1(T 5 
A/V 2 (white), respectively. The zero bias features correspond 
to Kondo peaks in dI/dV s <i- Black arrows indicate prominent 
inelastic features. When these inelastic features approach the 
boundaries of the Coulomb blockade region, the levels shift 
and alter lineshape (white arrows). Black dashed line in left 
map traces an inelastic feature across the Coulomb blockade 
region boundary and into the Kondo regime. (From Ref. 71) 



devices and semiconductor QDs. Shot noise is the result 
of current fluctuations that occur due to the discrete na- 
ture of charge carriers, and can provide a direct measure- 
ment of the effective charge (e*) in a system. For a com- 
pletely uncorrelated stochastic process, current is char- 
acterized by Poissonian shot noise. However, electron- 
electron interactions can correlate the charge transport 
process and suppress shot noise. 

In Djukic et ai, shot noise was measured in single 
molecule devices consisting of a D2 molecule between 
Pt source and drain electrodes^ Their measurements 
indicate that the shot noise is strongly suppressed with 
respect to the classical noise levels, and that transport 
through the Pt-D2-Pt system is carried predominantly 
by one nearly transparent conduction channel. Despite 
its Fermi liquid nature, many body effects, such as 
Coulomb repulsion and the Pauli exclusion principle, 
persist at low energies in the Kondo regime, and the 
effective charge is expected to be 5/3e due to these 
interactions 4 s - This is predicted to be a universal result 
in the Kondo system, independent of Tk- Measure- 
ments of the shot noise in the Kondo regime of a 
single molecule device would provide a test of this the- 
ory and contribute to the study of correlations induced 
in mesoscopic transport by different types of interactions. 

BEYOND THE SINGLE CHANNEL SPIN- 
\ SYSTEM 

We have heretofore considered primarily experimen- 
tal results concerning the single channel spin- 1/2 Kondo 
effect, as described by a localized Anderson impurity 
screened by delocalized conduction electrons. But the 
Kondo effect can occur when there is more than a single 
channel coupled to the impurity, and it can be observed 
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FIG. 16: a) Quantum phase transition depicting quantum state \A > driven by a non-thermal external parameter A to another 
quantum state \B > with a different symmetry, passing through a critical point at A = A c . In the SMT devices here, \A > is 
a singlet state and \B > is a triplet state that is partially screened by one conduction electron channel, represented by black 
arrows, b) Colormap over two Coulomb diamonds plotting dl/dV (in units of 2e 2 /h) as a function of bias voltage V and gate 
voltage V g at T — 35 mK and B — T. c) Details of the colormap in the dotted white rectangle in (c), showing the singlet to 
triplet spin transition. (Reproduced from Ref. 97 with permission from the Nature Publishing Company). 



due to the screening of a localized impurity with spin 
state S > 1/2. Moreover, it can result from non-spin de- 
grees of freedom, such as the SU(4) Kondo effect due to 
orbital degeneracies in carbon nanotube devices^ Until 
recently, experimental observations of these more exotic 
Kondo systems have been confined to interactions with 
semiconductor QDs and carbon nanotubes. 89 ! 90 This can 
be attributed, in part, to the limited ability with which 
one can tune the coupling strength between the active 
element and the source, drain, and gate electrodes in sin- 
gle molecule devices. The relatively high energy scales at 
which most SMT experiments have been conducted has 
also made the properties associated with these variant 
systems difficult to observe. 

This review has focused on the relatively common case 
of a single screening channel coupled to a single spin- 
1/2 impurity. As T OK this describes what may be 
called a fully screened Kondo effect. More generally, if 
the number of screening channels coupled to an impurity 
state is n and the impurity (molecule) is in spin state 5, 
the case of 25 = n corresponds to a completely screened 
Kondo effect. When n < 25 the impurity is said to 
be underscreened. Likewise, when n > 25 there are an 
excess number of screening channels and the impurity 
is overscreencd. The transport properties in the Kondo 
regime, namely the conductance measured as a function 
of temperature and bias, will depend on the spin state 
of the impurity and on the number of screening channels 
utilized by the conduction electrons. Recent reports have 
observed single molecule devices in even charge states 
exhibiting Kondo physics. 

Both Refs. 91 and 92 each involve devices containing 
a single Ceo molecule with 5 = 1, coupled to conduction 
electrons. This total spin arises because of the occupation 
of two different molecular orbitals each with one electron, 
with Hund's rule favoring the triplet ground state of to- 



tal spin 1. Each orbital can have different couplings to 
the conduction electrons in the source and drain, leading 
to two different energy scales, Tjci and Tft- 2 ' 93 ' 94 These 
energy scales depend exponentially on coupling strength, 
so that it is not difficult experimentally to be in the situ- 
ation where Tki <C Tki- If the experiments probe tem- 
peratures Tki €T< Tki, only one screening process is 
apparent, and the zero-bias resonance measured in Refs. 
91 and 92 corresponds to an underscreened spin-1 Kondo 
effect. 

The functional form of the equilibrium conductance 
as a function of temperature can be distinguished from 
the form used to describe the fully screened spin- 1/2 ef- 
fect by the more gradual slope and a value of G(T, 0) 
that does not saturate as T goes to K (Fig. I17h ). If 
temperatures below Tki can be reached then the second 
screening channel may allow for the impurity to become 
fully screened as T — > 0. A fully screened spin-1 impu- 
rity produces no resonance, so the curve of G(T, 0) will no 
longer increase monotonically with decreasing This 
can be described as a two-stage Kondo effect (Fig. 117b ). 

To date, few experiments with single molecule break 
junction devices have been carried out in the millikelvin 
temperature ranged where some of these non-spin-1/2 
Kondo effects may become more accessible. This is due 
in part to the large number of devices typically required 
for a successful analysis, and to the fact that the rel- 
atively high Kondo temperatures frequently associated 
with single molecule devices have made experiments in 
this temperature range an unnecessary complication. Re- 
cently however, in Roch et al^ investigators were able 
to conduct a detailed analysis of Kondo physics in Ceo- 
based electromigrated breakjunction devices by perform- 
ing transport measurements below 100 mK using a di- 
lution refrigerator. They observed multiple charge state 
transitions and the clear evolution from a Coulomb block- 
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FIG. 17: a) Equilibrium conductance vs. temperature mea- 
sured in the odd (red squares) and even (blue squares) charge 
states of sample A in Ref. 91. Solid and dashed lines are fits 
to the numerical renormalization group analysis for a fully 
screened S — 1/2 impurity and an underscreened S = 1 impu- 
rity, respectively. See reference for details, b) Two screening 
channels with different couplings to the local moment, and 
thus different Tk , can allow the conduction electrons to fully 
screen an S — 1 impurity for T < Tk i , leading to a zero-bias 
conductance that changes non-monotonically with tempera- 
ture. (Reproduced from Refs. 91 and 95 with permission 
from the American Physical Society). 

ade in even charge states to a Kondo resonance in odd 
charge states. The low energy regime in which their 
transport measurements are conducted enabled them to 
observe a number of novel effects. 

In one of their experiments, a detailed study of trans- 
port in an even charge state of a SMT reveals Kondo 
exchange interactions, evident in the stability diagram, 
attributed to higher spin (S > 1/2) Kondo states^ 
Most significantly, a quantum phase transition induced 
by gate voltage tuning is observed between two com- 
peting ground states (singlet and triplet electron spin 
states) involving Kondo correlations. This represents an 
observation distinctly different from previously observed 
Kondo transitions in semiconductor and nanotube quan- 
tum dots . 54 ' 98 In the singlet ground state no zero bias 
resonance is observed at base temperature. However, 
when the applied source-drain bias exceeds the energy 
difference between the singlet ground state and triplet 
(excited) state, a non-equilibrium Kondo effect is ob- 
served, producing dl / dV peaks at finite bias. This finite- 
bias Kondo effect has also been reported by Osorio et al. 
in transport measurements through a thiol end-capped 
oligophenylenevinylene molecule^ 9 - In the triplet ground 
state, a zero-bias peak due to Kondo screening is clearly 
visible, as are finite bias peaks related to the singlet 
(excited) state when the source-drain bias exceeds the 
triplet-singlet energy difference. 

A theory for a related gate-induced quantum phase 
transition in a SMT has been been put forth by Kirchner 
et al., in which a molecule between an antiparallel 
arrangement of ferromagnetic leads could couple to 
both the usual fermionic electron bath as well as the 



bosonic spin waves in the leadsi 1 ^ The authors propose 
that transport measurements in the Kondo regime 
could reveal the crossing of the critical point as a 
predicted change in the exponent of the fractional power 
law governing the behavior of the equilibrium Kondo 
conductance as a function of temperature. 

CONCLUSION 

Renewed interest in the highly correlated electron state 
leading to the Kondo effect has sparked a number of 
studies, both experimental and theoretical, addressing 
the manybody problem leading to this emergent phe- 
nomenon. First observed in an electrostatically defined 
GaAs QD, recent advances in nanofabrication techniques 
have enabled the experimental realization of the sin- 
gle impurity spin- 1/2 Kondo effect in molecular devices. 
Here we have reviewed experimental advances in the 
study of the Kondo state as observed in molecular elec- 
tronics. 

Transport measurements through individual molecules 
or molecular assemblies, using either a STM or a break 
junction configuration, have demonstrated the character- 
istic zero bias conductance resonance indicative of the 
Kondo effect. This resonance occurs due to an exchange 
process that results in the screening of the local mag- 
netic moment, represented by an unpaired spin on the 
molecule, by the sea of dclocalized conduction electrons. 
Single molecule devices represent a platform for prob- 
ing these interactions at more accessible temperatures 
compared to semiconductor QDs. Furthermore, the abil- 
ity to study transport in this regime using a variety of 
different systems {e.g. QDs, carbon nanotubes, single 
molecules) aids researchers in distinguishing between uni- 
versal behaviors and system related influences. Trans- 
port measurements in the Kondo regime of SMTs can 
exhibit characteristics not seen in their QD counterparts 
and not fully accounted for by commonly used theoretical 
models. This highlights the need for improved modeling 
techniques capable of incorporating complex correlations 
in a realistic molecular system. Additionally, it illus- 
trates the necessity for elucidating the unique qualities 
of the single molecule system, such as the precise way in 
which molecular orbital levels are renormalized when a 
molecule is coupled to metallic electrodes. Many open 
question remain in the study of Kondo physics in molec- 
ular electronic devices, leaving much room for future in- 
vestigations of this complex phenomenon. 
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